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ABSTRACT 

We present the results from new Nobeyama Milhmeter Array observations of 
CO(l-O), HCN(l-O), and 89-GHz continuum emissions toward NGC 604, known 
as the supergiant H ll region in a nearby galaxy M 33. Our high spatial resolution 
images (4'.'2 x 2'.'6, corresponding to 17 pc x 11 pc physical size) of CO emission 
allowed us to uncover ten individual molecular clouds that have masses of (0.8 - 
7.4) xIO^Mq and sizes of 5 - 29 pc, comparable to those of typical Galactic giant 
molecular clouds (GMCs). Moreover, we detected for the first time HCN emission 
in the two most massive clouds and 89 GHz continuum emission at the rims of 
the "Ha shells" . The HCN and 89 GHz continuum emission show offset from the 
CO peak and are distributed to the direction of the central cluster. Three out of 
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ten CO clouds are well correlated with the Ha shells both in spatial and velocity 
domains, implying an interaction between molecular gas and the expanding H ll 
region. The CO clouds show varieties in star formation efficiencies (SFEs), which 
are estimated from the 89-GHz and combination of Ha and Spitzer 24-/im data. 
Furthermore, we found that the SFEs decrease with increasing projected distance 
measured from the heart of the central OB star cluster in NGC 604, suggesting the 
radial changes in evolutionary stages of the molecular clouds in course of stellar 
cluster formation. Our results provide further supp ort to the pic t ure of sequential 



star formation in NGC604 initially proposed by iTosaki et al.l ( 120071 ) with the 
higher spatially resolved molecular clouds, in which an isotropic expansion of the 
H II region pushes gases outward and accumulates them to consecutively form 
dense molecular clouds, and then induces massive star formations. 

Subject headings: galaxies:individual (M 33) — ISM: individual objects: NGC 604 
— ISM: Hll regions — ISM: clouds — ISM: molecules — radio lines:ISM 



1. INTRODUCTION 



The expansion of an H ll region can trigger star formation in various ways (see the 
review of lElmegreenlll998l ) : direct compression of pre-existing clouds and globules due to high 
pressure of H ll regions ( "globule squeezing" ), and gas accumulation and cloud collapse at the 
edges of Hll regions and supernovae remnant ("collect and collapse"). The globule squeezing 
process works on the smaller scale of cometary globules, while the coUect-and-collapse process 
on the larger scale of the Hll region. Among above two processes, the collect and collapse 
process, first proposed by lElmegreen fc Ladal (Il977l ). is particularl y interesting as it allow s 
the formation of massive fragments to form massive stars or clusters (jWhitworth et al.lll994l ). 
As an Hll region expands with supersonic velocity, the ambient medium is photoionized and 
dynamically swept up into shocked layer, which may become gravitationally unstable and 
collapse to filaments, and then into one or more dense cores to form new generation stars. 
Many multi-wavelength studies for the collect and collaps e process have b een performed 
in our Galaxy and neighboring galaxies (see examples in lElme green! Il998l ). From JHK 
photometry. The stellar clusters of the different age are seen spatially resolved in an H ll 
region, and more over, embedded younger massive clusters are found to lie at the perip heries of 



Hll regions (e.g.. iDeharveng et al.ll2005l : IZavagno et al.ll2006l : iDeharveng et al.ll2008l ). These 



clusters are possible second-generation clusters whose f ormation has be en triggered by the 
expansion of Hll region ("propagating star formation"; lElmegreeru 119981 ) . 
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Mayyal ( 119941 . Il995[ ) and iMayya fc Prabhul (119961 ) extended this process of propagating 
star formation in molecular clouds to a larger scale to explain the star formation process in 
Giant Hll regions (GHRs). GHRs can easily spread over areas 10^-10^ pc across, and play 
a vital role in determining the large-scale structures of their host galaxies. However, little 
is known about the possible connection between ambient interstellar matter and expansion 
of GHRs over such distances. This is because there is not such a giant Hll region in our 
Galaxy, and thus toward external galaxies it has been difficult to perform high-resolution, 
high-sensitivity, and wide-field observations at once. Evidences of "collect and collapse" 
process on large-scale were r nost notably frorn the correspondence between th e old H I holes 



and young OB associations (IBrinks fc Bajajalll986l : iDeul fc den 



iartoei 119901). showi i ig the 



kinematic and spatial evidences of expansion of H I and Ha shells. iDeul fc den Hartogl (Il990[ ) 
allow for the possibility that many such holes have been excavated by massive star formation 
activity and formed by the local collapse of H I atomic gas to form giant molecular clouds 
(GMCs) which are the principal sites of star formation . In fact, the GMCs tend to be 
spatially clustered around the Hi holes (lEngargiola et al.ll2003l ). 



The spatially resolved studies of these GMCs are necessary in order to understand the 
triggers of star formation around GHRs through the formation of dense molecular clouds 
which will be taken place in individual GMC. Recent progress of millimeter arrays provides 
wide-field, high-sensitivity, and high-resolution images of the molecular gas for relatively 
shorter observational time, which has potential to reveal clumps within GMCs and physical 
relation to ongoing star formation. Our objective is to study how extragalactic GMCs evolve 
to form stars through dense gas formation which might be affected by large-scale phenomena 
that is, expansion of giant Hll region. In particular we focus on the evolution of GMCs by 
probing their physical nature such as their density structures and star formation activities, 
and their spatial distributions against the activities in galactic environments from the multi- 
wavelength observations. With this purpose, we embarked upon detailed studies toward 
nearby extragalactic GMCs in a bright supergiant H ll region, NGC 604. 

The giant Hll region NGC 604 is located in northern arm in the spiral galaxy M 33. 
The prox imity (distance of 0.84 M pc; iFreedman et al.ll200ll ) and favorable inclination of M33 
{i = 51°: ICorbelli fc Saluccill2000l ) are sufficient to resolve gas components with little contam- 
ination along the line of sight and to under stand the inner structure of GMCs (i.e., typically 
40 pc in diameter; IScoville fc Sanderslll987l ). Because of its striking fea tures, the giant H ll re- 



Wilson fc Scoville 11992 



gion NGC 604 has been observed in a wide range o f wavelengths; radio (jViallefond et al.l 



Churchwell fc Goss 



1999 



Tosaki et al.ll2007l). optical ( 



1992 



Drissen et al. 



19931: 



Hunter et al.l 



Tiillmann et al 



1996 



Tenorio-Tagle et al. 



2000 



Maiz-Apellaniz et al.l 120041 ). and X-rays 



200^ 



sl). According to previous stu dies, the Ha ne bula has a core-halo struc- 
ture extending up to a radius of 200 - 400 pc (iMelnickl Il980l . see also Figure |6]). The 
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Ha core, containing a very young (3-5 Myr) and rich (over 200) 0- type stellar popula- 
tion, is surrounded by photo-ionized filaments and a section of shells 



Diaz et al 



2000 



1996 



Hunter et al 



Bruhweiler et al 



2003 



1996 



Gonzalez Delgado fc PerezI 12000 



Drissen et al. 


1993; 


Tenorio-Taele et al. 



Relano fc Kennicuttll2009[ ). The energy of the stellar wind has 



been estimated to be enough to form the core of the Ha nebula (e.g., iTenorio-Tagle et al. 
2000l : iMaiz-Apellaniz et al.ll2004j : iTiillmann et al.ll2008l ). on the other hand, SNRs as we ll as 
stellar winds would contribute to blow it up to the halo structure (e.g.. lYang et al.lll996l ). As 
noted in Figure El we call the shell surrounding the main cavity "Shell A" and the shell lo- 
cated to the south of Shell A, "Shell B" , which a ppears to be a blowout through the molecular 
tunnel over Shell A (ITenorio-Tagle et al.ll2000[ ). OB and W-R stars are mostly distribute d 
in Shell A, while, a small number of OB stars in Shell B (see Figure 7 in lYang et al.l 119961). 
Near-i nfrared JHK photometric study of NGC 604 was recently reported in iFarina et al. 
( 120101 ). the candidates of massive young stellar objects (MYSO) appear aligned with bright 
Ha Shells A and B. An Hi hole surrounds around diffuse Ha nebula and some GMCs are 



within 150 pc around the Hi hole edge (lEngargiola et al.ll2003t iHeiner et al.l 120091 ). The ex 



istence of five GMCs over t^ 
previous ^^CO observations ( 


le extent Ha emission in the NGC 604 region was revealed by 


Viallefond et al. 


1992; 


Wilson & Scoville 


1992; 


Engargiola et al. 


2003; 


Rosolowskv et al. 


200c 


1). The sensitivity and spatial resolution of these observations 



using the OVRO interferometer and the BIMA array were, however, not sufficient to resolve 
the inner structure of any of these CMC. 

We present here the new results of the interferometric CO(J=1-0), HCN, and 89 GHz 
continuum observations toward the NGC 604 region. Our imaging achieved a higher spatial 
resolution than any other studies so far (10 pc scale for ^^CO emission). We observed the 
J=l-0 emission lines of ^^CO and HCN, and adjacent 89 GHz continuum emission. These 
three probes trace different components of the CMC; i.e., relatively diffuse molecular clouds 
(the critical density of ^^CO line emission, ~ 10^ cm~^), denser regions (riHa > 10^ cm~^), 
and active star-forming regions, respectively. Indeed, the HCN line emission has been known 
to trace de nser regions in mole cular clouds linked to star- formation activities rather than 



12 



CO (e.g., ISolomon et al.lll992[ ). The 89 GHz continuum is considered to be mainly domi- 
nated with free-free emission from ionized gas by young OB stars, and thus expected to trace 
massive star forming region. The HCN observation and millimeter-wavelength continuum 
emission toward NGC 604 are reported for the first time in this work. Besides, this work 
provides the first example of extragalactic GMCs in spiral galaxies being observed on a size 
scale comparable to typical galactic GMCs. 

The outline of this paper is as follows. The observations and data reduction are summa- 
rized in §|2j In §[3] the general procedure for identifying clouds in the survey and the results 
of each probe are described. In the discussion in § HI we quantify the star formation activities 
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in the GMCs using our data and also 24 /im and Ha data, and discuss the varieties of the 
physical states of the identified clouds. We also discuss the relationship between the nature 
of these clouds and their evolutionary stage. Finally, we interpret our results i n the context 



of a s equential star formation scenario, as initially proposed for NGC604 by iTosaki et al. 



(120071 ) from the spatial relationship between the GMCs and Ha nebula. The summary of 
this work is given in ^ 



2. Observations and Data Reduction 

We performed ^^CO(l-O), HCN(l-O), and 89 GHz continuum observations toward NGC 604 
using the Nobeyama Millimeter Array (NMA) at Nobeyama Radio Observatory (NRO) from 
November 2004 to December 2006. We observed ^^CO line emission toward the five fields of 
view (FoVs) with both C and D configurations of NMA (see Figure [H Tabled and Table [2]). 
The area of the five FoVs was about 110" x 145", corresponding to 450pcx600pc covering 
the entire spatial extent of the optical nebulae. Besides, HCN(l-O) and 89 GHz continuum 
emission were observed toward two out of the five FoVs to cover the two largest GMCs. To- 
ward other three FoVs including smaller GMCs, we have not observed because it is expected 
to need more integration time to detect than toward the two largest GMCs. The sizes of 
FoVs are represented by the full width at half-maximum (FWHM) of the primary beam at 
each frequency and they are given in Table [H 

We used DSB SIS receivers as the front-end, and the digital spectro-correlator New-FX 
(1024 channels, 32 MHz bandwidth) as the back-end for the ^^CO and HCN observations. 
The total bandwidth of the FX correlator corresponds to 83.2 km s~^ at the ^^CO frequency, 
and 121.6 km s~^ at HCN, respectively. The 89 GHz continuum observations were performed 
simultaneously with the HCN observations using the Ultra- Wide-Band Correlator (UWBC; 



Okumura et al.ll2000l ) which has 128 channels and a bandwidth of 1024 MHz. The continuum 
emission data were made by merging the frequency information in the IF signals from lower 
sideband (LSB) of the UWBC. The frequency ranges containing line emissions of HCN 
and HCO^ (89.188526 GHz) were removed from the original visibility data. The effective 
bandwidth of line-free channels is 768 MHz. 

The bandpass calibration was done using the quasar 3C 454.3. The radio sources 
B0234+285 and B0133+476 were observed every 20 minutes as phase calibrators. The flux 
scales of these calibrators were determined by comparison with Uranus and Neptune. The 
uncertainty in absolute flux scale is expected to be less than 15%. 



The interferometric data were reduced using the NRO software package UVPROC-II 



- 6 - 



( iTsutsumi et al.lll997l ). The final mosaiced image was made by combining the visibihty data 
of all the available FoVs with the software package MIRIAD (using the task invert to make 
the synthesis image, and mossdi to CLEAN the image for the line data or mfclean for the 
continuum data). All these maps were made with natural uv weighting. The synthesized 
beam sizes of the ^^CO images for each FoV are given in Table [H while those of the HCN 
and 89 GHz continuum are in Table [2l Note that the spatial resolutions in the partially 
overlapped areas between FoVs are the average of those in the overlapping FoVs. 

The final 3-dimensional data cube of the ^^CO and the HCN data has a frequency reso- 
lution of 1 MHz (corresponding to 2.6kms~^ at 115 GHz, 3.4kms~^ at 89 GHz, respectively) 
and the typical noise level of them are given in Table [1] and Table [2j The noise level of the 
89 GHz continuum data is also given in Table |5J The noise level was measured over the 
emission-free area. The final flux densities and noise levels for the sources were estimated 
after correcting for the primary beam attenuation. Note that the noise levels rise by a factor 
of 1.5 at the edges of each FoV due to the primary beam correction. Since the uv coverage 
of each pointing is not identical, the synthesized beam patterns differ among pointings. We 
briefly summarize the observational parameters in Table[T]and Table[2]for the interferometric 
observations. 

It is worth noting that interferometric observations do not have the sensitivity to large 
scale features (i.e., low spatial frequencies) of size larger than ~ A/i?min, where A is the 
observing wavelength and -Bmin is the minimum projected baseline. For our observations, 
-Bmin is 8-9 kA at 115 GHz, and 4 kA at 89 GHz. Thus emission structures extended over 
spatial scales more than 23"-26" (93-105 pc) in 115 GHz observation, and 52"(~ 211 pc) in 
89 GHz observation are resolved out in the interferometric images. In order to estimate the 
missing flux in the CO observations, we compared the total integrated fluxes for each FoV 
in our NMA maps to the single-dish maps obtained using the 45-m telescope at NRO (R. 
Miura et al., 2010, in preparation). The fluxes were measured within a box of 42" x 42" 
at the center of each FoV. The interferometric to single-dish flux ratios are given in column 
8 of Table [1] The estimated ratios of FoV2 and FoV4 are relatively lower than other three 
FoVs, less than 30%, because there is little molecular gas that detected in interferometric CO 
observation although the single-dish maps in FoV2 and FoV4 includes the envelope of the 
three largest GMCs which are located at the center of FoVl, FoV3, and FoV5, respectively. 
The averaged ratio over the five FoVs is 0.46 ± 0.17, suggesting that half of the CO emission 
is resolved out in our interferometer observations. 
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Results 



3.1. i2co(l-0) 

The CO integrated intensity map in the velocity range Vlsr = —211 km to —256 km 
is shown in Figure |2j We identified CO clouds using the following criteria; cloud emission 
must be more than three times stronger than the rms noise level (cxrms) and at least two 
times wider in velocity than our velocity resolution, and the cloud sizes must be larger than 
the equivalent area of the synthesized beam size. Ten clouds were identified with these above 
criteria. Their physical properties were obtained as described below. 

The cloud size D is defined as the average of the deconvolved FWHM diameters in 
right ascension and declination. The typical error of D was estimated as an average of the 
spatial resolution in the individual pointings. The velocity FWHM (Vfwhm) was derived 
from Gaussian fits to emission line. Virial masses ( Myjr) and molecula r cloud masses (Mmoi) 
were calculated using the procedures described in I Wilson &: Scovilld ( Il990l . and references 
therein). They were derived using the equations, Mvir = 99DpcV|^jjj^ and Mmoi = 
1.61 X 10^(i^pj,S'co Mq, where D^^ is the cloud size in parsec, c^Mpc is the distance to M 33 
in Mpc, and Sqq is the CO fiux density in Jykms"^. We used 1.4 times the deconvolved 
diameter of the cloud in calculating the viria l mass, since this diamet er was found empirically 
to contain most of the fiux of the cloud (jWilson &: Scovilld Il990l ). Here we assume the 
galactic value of the conversion factor from CO fiux densities and H2 column densities, 
(3 ± 1) X 10^''cm~^/(Kkms~^) ( IStrong et al.l 119881 ). and a factor of 1.36 has been included 
for the contribution of helium. 

Table [3] summarizes the properties of the ten clouds we identified. The positions of the 
ten clouds are shown in Figure [2l The cloud sizes range from less than 5 pc to 29 pc and their 
masses are in the range ( 0.8 - 7.4) x 10^ Mfr^, which are comparable to those of typical GMCs 



in the Milky Way (e.g., IScoville fc SandersI 119871 ). Six clouds (NMA-1, NMA-3, NMA-4 



NMA-5, NMA-6, and NMA-10 in Table [3]) were ne wly identified with our observations. The 



rest w ere already reported by the previous papers (IWilson fc Scovilldll992uEngargiola et al. 
2003h . 



We would like to mention two other clouds, NMA-11 and NMA-12. Even though their 
features do not satisfy our criteria for cloud identification, they were detected over 3 a but 
only one velocity channel, they were also detected in our single dish image with the NRO 
45-m telescope (R. Miura et al., 2010, in preparation). One source NMA-11 is located to the 
north of NMA-8, while the other source NMA-12 is to the north of NMA-9. The properties 
of these two clouds are also given in Table [31 
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3.2. HCN 

Two HCN emission components were marginally detected more than three times stronger 
than the rms noise level and over two velocity channels (each has 3.4 kms~^). One is at 
(01^34°^33!4, +30°46'50'.'5) at the velocity of -245 km s'^ and the other at (01^34'^34!3, 
+30°46'21'.'0) at the velocity of — 225kms~^. Figure H] shows the HCN integrated intensity 
map overlaid on the CO map, which were integrated over the velocity range where the emis- 
sion was seen (> 3 arms)- These two HCN clouds seem to be associated with CO clouds, 
NMA-8 and NMA-9. The peaks of both HCN clouds show offset from those of CO clouds, 
and moreover, both of HCN clouds are located on the northern side of each CO cloud, to 
the direction of center of H ll region. To check this association we have compared the veloc- 
ities of the CO clouds to those of the HCN clouds. The HCN and CO spectra of these two 
components at the velocity resolution of 3.4 kms~^ are shown in Figured Each spectrum 
was measured over the area where the emission was seen (> 3o"rms)- The CO images were 
convolved to the same spatial and velocity resolution of the HCN images and the CO spectra 
were measured in the same region as the HCN spectra. The peak of the HCN spectra is 
slightly blue-shifted compared to that of the CO spectra in the same spatial region, which are 
found in both profiles of NMA-8 and NMA-9. The bluer velocity corresponds to the direction 
of downstream by considering the direction of galactic circular rotation (from south to north 
in Figure in this region, and hence these two HCN components seem to be associated with 
the blue components of CO clouds, which are located on the downstream (northern) side. 

This suggests that the components with different velocities have differe nt flux ratio of 



HCN t o CO. The flux ratio of HCN to CO is often used as a density probe (e.g.. lHelfer fc Blitz 



1997al ). We convolved the CO images with a Gaussian function of the same beam size as 
that of the HCN images for comparison between the total flux densities in the CO and HCN 
images at the same spatial scale. We obtained each total flux density by integrating over 
the extent of the emission components with more than 3o"rms- The flux ratios of HCN to 
CO (Shcn/5£o 0) were estimated to be 0.024 ± 0.007 for NMA-8 and 0.018 ± 0.007 for 
NMA-9. For the CO clouds where we did not detect any HCN emission (NMA-1, NMA-4, 
NMA-6, NMA-7, and NMA-10), we set upper limits on the line ratio using the detection 
limit of 3crrms- The upper limits of Sylc^^/Sqq ratios are given in Table HI 



There is similar CMC's HCN survey in M 31 using IRAM 30-m radiotelescope (iBrouillet et al 



20051 ). We found that our estimates of Sucn/Sqq are comparable @to ratios of 0.0075-0.028 



^We define the CO flux densities in tlie convolved images as Sqq. 

^This might infer that the missing fluxes in our interferometric observations are nearly equal between the 
CO images and the HCN images 
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that obtained from their spectrum data of M 31 at a spatial scale of ~100pc (~28"). The 
S-RCN / Sqq ratios of the two HCN coi nponents are also corr iparable to those of typical galac- 
tic GMC (0.014 ± 0.020 on average; iHelfer fc Blitall997bl ). which were measured in nearby 
GMCs at the spatial resolution of 0.5 pc. 



3.3. 89 GHz continuum 

We detected the 89 GHz continuum emission in NGC 604 for the first time. Four con- 
tinuum components were detected above 4-crrms noise level in our synthesized image. The 
two strongest detections (> 5 Urms) in our NMA images are labeled as Source I and Source II. 
The right panel of Figure [5] shows the 89 GHz continuum map overlaid on the CO integrated- 
intensity map. 

Source I was detected over 6 arms at the position (0l''34™'33':5, -|-30°46'55'.'9) and appears 
to be associated with NMA-8 (within 9"). Source II was detected over 5 cirms at the position 
(01^34'°32n, +30°46'59'.'3) which is consistent with the position of the CO emission peak of 
NMA-4. Thus we consider that these two 89 GHz components can be physically associated 
with NMA-8 and NMA-4, respectively. Moreover, we found that the distribution of these 
sources. Source I and Source II, exhibits a similar i ty to the 8.4 GHz continuum emission 



components A-D presented in IChurchwell Sz GossI fll999l ): our Source I corresponds to the 



position of their A and B components, while our Source II to the position of their C and D 
components (see Figure [5]). The other two components. Source III and Source IV, are located 
at (01^34"^32!4, +30M7'18'.'8) and (01^34^32^8, +30°46'28'.'8), respectively. The signal-to- 
noise ratio of these components is low (< ScTrms) and Source III is located at the edge of the 
FoV. In addition, neither CO clouds nor 8.4 GHz emissions are found to be associated with 
them. Therefore, we do not consider the Sources III and Source IV for further discussions. 

We measured the flux densities over the region where the emission was seen (> 3 arms) 
and the errors of the total flux densities were estimated taking account the uncertainties 
in the absolute flux scale and the noise fluctuations. The estimated flux densities of the 
two components are 3.4 ± 0.9 mJy and 9.1±1.5mJy in the NMA-4 region (Source II) 
and in the NMA-8 region (Source I), respectively. The millimeter- wavelength continuum 
flux provides a constraint on the emission measure of ionized gas, from which we can ob- 
tain information about the number of stars present for models of number of UV photons. 
This is because the free-free emission from the gas ionized by OB stars contributes largely 
to millimeter- wavelengths and is optically thin at millimeter- wavelengths. Other contribu- 
tions to millimeter-wavelength continuum, non-thermal radiation from supernova remnants 
(SNRs) and thermal emission from interstellar dust, are expected to be small. Assuming 
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that the 89 GHz continuum flux densities are entirely dominated by free-free emission from 
the ionized gas, we calculated the total production rate of ionizing photons from the observed 
89 GHz continuum flux densities. 



In the calculation, we followed equation (2) in lScoville et al.l ( 1l99ll ). which provides the 
relationship of the free-free emission flux density at 110 GHz to the total production rate of 
ionizing photor is for optica l ly th in plasma. Using the thermal (free-free) spectral index of 
—0.1 in M 33 (IBuczilowskil Il988l ) . the relationship of the free-free emission flux density at 
89 GHz (S'ggGHz) and the total production rate of ionizing photons (Q*s~^) is given by: 



S> 



89GHz 



^.58 X 10" 



,Q*/lO 



49^,-1 



"Mpc 



-mJy, 



where c^Mpc is the distance to M 33. The Hll regions are assumed to be ionization-bounded 
at a temperature of 10^ K. The total production rate of ionizing photons are (28.0 ± 7.3) x 
lO'^^s"^ in the NMA-4 region (Source 11) and (74.8 ± 12.3) x 10^^ s"^ in the NMA-8 region 
(Source I). The value for Source I is consisten t with the total value expe cted from the 8.4 GHz 



fluxes of components A and B (68x10^^ s ^; IChurchwell fc Gosslll999l ). As for Source II it is 



not consistent with C and D components (64x10^^ s~^). It is possible that more flux from 
large scale structures in NMA-4 are resolved out than those in NMA-8, because missing flux 
in interferometric observations depends on the sizes of the emission sources. Assuming that 
the missing flux at 89 GHz continuum is the same as that of the CO images, the ionizing 
photon production rate for the NMA-4 region is consistent with that from the 8.4 GHz fluxes 
of components C and D. 



3.4. Overall features and Comparison with Ha 



In this section we compare the distributions of CO, HCN, and 89 GHz continuum emis- 
sions with the morphology of the Ha shells in NGC 604. The composite image of CO and 
HCN line emissions, and 89 GHz continuum emission overlaid onto the Ha image obtained 
by the Hubble Space Telescope (HST) is shown in Figure [61 

Three out of ten identified clouds (and NMA-11) are associated with the bright shells, 
i.e.. Shell A and Shell B (see Figure [6]). The bright western part of Shell A shows an 
hourglass-like structure, which appears to be fitted by the distribution o f ^^CO emission o f 
NMA-4. The velocity structures of these Ha shells show true expansion (lYang et al.lll996l ). 
possibly powered by stellar winds from the central OB star clusters and/or SNRs (Section 1). 
The line-of-sight velocity of NMA-4 (— 247kms~ ^), which correspond s to th at of the Ha shell 
at the same position (Ha emission reported in iTenorio-Tagle et al.l (120001 )). corresponding 
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to Shell A. Therefore, it is considered that NMA-4 is a cloud component associated with 
the expanding Ha shell, Shell A. The cloud NMA-11, the marginally detected cloud, is also 
located on the eastern part of Shell A. 

The cloud NMA-8 itself has an arc-like structure and is distributed to surround the 
cavity, nicely tracing the bright rim of Shell B. The 89 GHz continuum and HCN emission, 
which are considered to be associated with NMA-8, were detected in the inner wall of Shell B, 
i.e. to the direction of the central stellar cluster. Moreover, there is an offset between peak 
positions of three probes: CO, HCN, and 89 GHz continuum, which are aligned to the 

it velocity of NMA- 



direction of central stellar clusters (from south to north). The line-of-sig! 
8 (— 243kms~^) corresponds to that of the Ha shell at this position (ITenorio-Tagle et al. 



2000l ). corresponding to Shell B. Therefore, it is considered that NMA-8 is a cloud component 



associated with the expanding shell B. 

Looking toward the halo structure of the Ha nebula, the distribution of ^^CO emission 
of NMA-2, NMA-6, NMA-9, and NMA-12 do not overlap with the area seen in Ha emission 
but extend along the edge of diffuse Ha emission. Although it is difficult to address whether 
or not they are physically associated with the diffuse Ha emission due to lack of the spectral 
data of Ha, we would like to mention the interesting feature of the dense-gas-forming region: 
the HCN component detected around the NMA-9 cloud, showing offset between CO peak 
and HCN peak, is located in the northern part of NMA-9; i.e., to the direction of the Ha 
nebula. As for the rest of the clouds, they are seen on the diffuse Ha emission. The cloud 
NMA-10 is located in the southern part of the " U" -shaped Ha loop, which marginally shows 



expansion signatures (ITenorio-Tagle et al.ll2000l ). 



Consequently, we found the distribution of CO, HCN, and 89 GHz continuum emissions 
which we detected in our observations are well correlated with the expanding Ha shells 
in NGC 604, in both morphology and kinematics. In contrast with the concentration of 
molecular gas on the rim of shells, no molecular clouds are observed inside Shells A and B. 



4. Discussion 



4.1. Comparison with 24 /im and Ha emission 



In this section, we derive the star formation efficiencies (SFEs) from the available Spitzer 
Archival 24 /xm data and Ha data in order to discuss the star formation activities over all 
five FoVs, because a combi nation of Ha and 24 nm luminosities is a better tracer of the total 
star formation rate (SFR, iKennicutt et al.l 120071 ). We then compared the SFEs with those 
in inferred from our 89 GHz data. 



-VI- 



The 24 /iin emission is often used as an indicator of the current star formation because it 



is rad iated from dust grains heated by newly formed massive stars (e.g.. lRelafio fc Kennicutt 
20091 ). The Ha emission is also directly linked to star formation, however, it is only sensitive 
to the ionizing photons unabsorbed by dust. There is empirically known to be a linear 
combination of the 24 /zm emission, tracing the obscured star formation, and the observed 
Ha luminosity, tracing the un-absorbed star formation, which correlates bett er than other 



SFR tracers with the extinction-corrected Ha luminosity (jCalzetti et al.l 120071 ). 



We retrieved the 24 /zm images of M 33 obt ained with MIPS (IRieke et al.l l2004j ) from 
the Spitzer Space Telescope fjWerner et al.ll2004f ) data archive. The Basic Calibrated Data 
(BCD) was created using the Spitzer Science Center (SSC) pipeline. A background matching 
was applied between overlapping fields of view and background gradient was removed. The 
FWHM of the point spread function in the image is 5'.'7 and the grid size of the final image 
is 2'.'5 per pixel. Our 89 GHz continuum image was overlaid on the resultant 24 /im image as 
shown in Figure [51 Both Source I and Source II in our 89 GHz continuum image correspond 
to peaks in the 24 fim map. This is one of the evidences that the 24 /im and the 89 GHz 
continuum emission both trace sites of star formation embedded in dust clouds. 

The reduced data of Ha emission in M 33 were kindly provided by R. Walterbos. It was 
observed with the 0.6-m Burrell-Schmidt telescope at Kitt Peak National Observatory. The 
dimensions of the CCD are 2048x2048 with pixels of 2'.'03 and a total field of view of about 
70' X 70', and then trimmed with t he same size of our C O image in NGC 604. The sensitivity 
is 0.8 X lO^^ergs"^ cm~^ arcsec"^ (iHoopes et al. II2OO1I). More details of the observations and 
reduction process are described in iHoopes &: Walterbos (2000). 



4-.1.1. SFE from 89 GHz continuum emission 

The total stellar mass can be compared with the mass of the molecular clouds to measure 
SFE. The SFEggcHz is defined as follows: 

'Si* i^SQGHz - TT ~T7^' ' ' 

where M^tar = Jj^j.^''N{M)MdM is total stellar mass between Mmin and Mmax, and N{M) 
represents stellar mass function. We defined the missing-flux-corrected molecular cloud mass 
as Mj^qj using the missing flux derived in Section [2l We analyzed total number of stars from 
the total ionizing photon rates (Secti on 13. 3p using the ionizing rates of each 0-type stars 



(09.5 V - 03 V iMartins et al.ll2005[ ). because these stellar spectral types are considered 



to be largely responsible for the ionizing photons of interstellar matter. We assumed that 
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the Hll region NGC 604 have Salpeter IMFs (2.35; iGonzalez Delgado fc PerezI 120001) with 
mass hmits of 0.1-58 Mq (corresponds up to spectral type 03 V; iMartins et aLll2005l ). The 
number of 0-type stars (09.5 V - 03 V) in the NMA-4 and NMA-8 region is 185±31 and 
70 ± 18 stars, respectively. The SFEggcHz was estimated to be 12.4 ± 6.6 % for NMA-4 and 
5.1 ± 2.3 % for NMA-8. As for other clouds which were not detected in continuum emission, 
upper limits of SFEggcHz were calculated and they are given in Table HI 



4-1-2. SFE from 24lJ'rn emission and Ha emission 



The Ha and 24 /im images were convolved into a common 5'.' 7 resolution and regridded 
to 2'.'5 per pixel. We measured the 24 /xm and Ha luminosity over the same spatial regions 
of all the CO clouds we identified. We measured the integrated flux over the extent of H a 
nebula and comp ared with previously reporte d fluxes from the literature (jKennicuttl Il984l ) . 
The Ha fluxes in iRelano fc KennicuttI toO^ are computed fro m the observed Ha +[N II] 
fluxes using the spectroscopically determined [N II] /Ha ratios (IBosch et al.l |2002| ). The 
extinctioncorrected Ha luminosity is then expressed as a linear combinatioi i of these two 
luminosities: Lcorr(Ha) = Lobs(Ha) + a x L(24/im), where a = (0.031 ±0.006) fICalzetti et al. 



20071 ). Ionizing photon rates are derived from extinction-corrected Ha luminosities L(Hac 
using 

Q*{s~') = 7.32 X 10"L,o„(Ha) (ergs'^). 



(3) 



which assumes an ionization bounded case B nebula (IMayya fc Prabhulll996l . and reference 
therein). The star formation efficiency from 24 /xm and Ha (SFEHo,24/xm) was calculated 
using the same procedures described in previous subsection 4.1.1. The SFEHa,24Ajm for each 
cloud is given in Table HI The NMA-4 has the highest SFE of all clouds and there is the 
same tendency in both SFEs. The derived SFEHa,24^(m of 0.2 - 11.5% is good agreement 
with that from 89 GHz continuum (< 1.2 - 12.4%) within errors. 

The total amount of ionizing photon rate {Q*) in the NGC 604 region as a whole (within 
a radius of 400 pc) is 413 photons s~^. The summation of the Q* over identified clouds is 
188 xlO^^ photons s~^, which occupies about 46% of the total Q* of a whole region. On 
the other hand, the Q* from the central stellar clusters (~ 200 O stars) is estimated to be 
145 X 10^^ photons s~^, assurn ing that the central stellar clusters have the Salpeter IMF 
(2.35) and lMartins et al.l (120051 ) 's ionizing photon rates for each stellar spectral type, which 
contributes about 35 % to the total Q* of a whole region. Thus this implies that such a large 
ionized region has been created not only because of the central stellar cluster but also the 
embedded stellar clusters in the clouds. 
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4-1.3. Comparison to other H II regions 



In comparison, we derive SFEs of three well-known Hll regions (Orion A, W49A, and 30 
Doradus) in the Local Group. The first two Hll regions contain similar molecular cloud mass 
(i.e., mass of l-TxlO^M©) to our identified CO clouds in NGC604. Orion A is the nearest 
and most well-studied region of star formation in our Galaxy, and is known to be excited 
by a few late 0-type star or early B-type star (e.g., ICardelli fc Clayton! 119881 ). W49A is 
one of the most massive and luminous H ll regions in our Galaxy and contain s 30 individual 
ultracompact H ll regions, corresponding to about that number of OB stars ( IConti fc Blum 
2OO2I ). The most luminous in the Large Magellanic Cloud (LMC), 30 Doradus giant Hll region 
is introduced to compare with the secondary luminous Hll region, NGC 604. Very luminous 
NGC 604-class H ll regions is absent in our Galaxy, and hence it is frequently compared to 
the NGC 604, although 30 Doradus has about 3-4 times as total ionizing flux as NGC 604. 
The estimated SFEs with the same assumption ab ove are tabulated in Table |5l The Q* oi 
Orion A was determined from the stellar types in ICardelli fc Clayton! (jl988[ ) and the SEE 
of Orion A is 2.2 %. The Q* of W49A was calculated from the measure ment of conabiniii g 
the 25/im luminosity (jWard-Thompson fc Robsonlll990l ) and Ha emission (jKennicuttlll984l ). 
and this can be a v alue o f 182x10^^ s~^. The Q* of 30 Doraus is derived only from Ha 



emission 



KennicuttI ( 1l984l ). and thus the SEE of 30 Doradus can be more than 52.8%. The 



SEE of NMA-4 is much larger than Orion A and comparable to those of the galactic Hll 
region W49A, while the SEE of NMA-8 is comparable to or slightly less than W49A. The 
SEE of each decomposed cloud in supergiant H ll region NGC 604 is quite similar to that of 
galactic H ll region. Both of NMA-4 and NMA-8 have about three to ten times lower SEEs 
than 30 Doradus, although they are needed to be compared on the same scale. 

The averaged efficiency of star formation out of giant molecular associations (GMAs; 



e.g.. lRand fc Kulkarnilll990l ) in other extragal actic GHRs tends to be a bout 1-10 % out of the 
parental molecular cloud mass (lOMO^MO, iMavva fc Prabhulll996r ). In order to compare 
SEEs of them with that of NGC 604 on the same scale, the total (missing-fiux-corrected) 
mass in NGC 604 region is estimated to be 51 xIO^Mq, and the extinction-corrected Ha 
luminosities Lcorr(Ha) is to be 4.3 x 10^^ s~^. The derived SEE over NGC 604 region is 7%. 
This is good agreement of other extragalactic GHRs. On the other hand, the SEE of the 
most luminous H II region 30 Doradus is extremely higher than NGC 604 or any other GHRs, 
which suggests its half of parental molecular clouds have been converted into stars and also 
dozens of 0-stars would have dissipated a whole GMC. NGC 604 may have potential to 
continue star formation activity, which leads to a 30 Doradus-scale GHR, since NGC604 has 
still contains a considerable reservoir of molecular gas to be used for star formation. 
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4.2. The Different Properties and Evolution of GMCs 

As we mentioned in the former section, we found that the identified clouds have different 
SFEs (see Table H). Especially, NMA-4, NMA-8, and NMA-9 have conspicuous differences 
in their SFEHQ-,24Atm- In this section, we focus on these clouds and discuss the relationship 
between their different properties and their evolutionary stages. 



4.2.1. NMA-4 and NMA-8 

NMA-4 is located near the central stellar clusters (the projected distance from the 
clusters is ~32pc) and also along the bright rim of Shell A (western part of Shell A; see 
also Figure E]) . The cavity surrounded by Shell A is known t o be filled with hot coronal 



gas that emits X-ray s and is transparent to the ionizing UV (jMafz-Apellaniz et al.l 12004 



Tiillmann et al.ll2008l ). The 8 fim emission, which is often used as tracer of photo-dissociation 



regio n, is mostly distributed around the UV emission (see Figure 2 in iRelano &: Kennicutt 



20091 ) ■ where the clouds NMA-8 and NMA-4 are located. The NMA-4 and a part of NMA-8 



are considered to be directly exposed to the UV radiation from young massive stars. 

In comparing the physical properties of NMA-4 with NMA-8, we found significant 
enough to distinguish between their evolutionary stages. The clouds NMA-4 and NMA- 
8 have the highest SEE and the second-highest SFE, respectively, and nevertheless, they 
are different by a factor of ~ 3. Its significant difference in SFEs is mainly due to their 
respective masses of molecular gas; NMA-4 has the cloud mass of ~0.8xl0^ Mq and size 
of ~10pc, while NMA-8 has the cloud mass of ~7.4xlO^M0 and size of ~29pc. These 
results might suggest that most of the cloud mass in NMA-4 has been more affected with 
star formation activities, i.e., conversion into stars, dissipation by the stellar winds, or dis- 
ruption by intensive photo-dissociation, than in NMA-8. The cloud NMA-8 is one of the 
largest molecular clouds and shows the obvious sign of dense gas formation (Section 13. 2p . 
Since NMA-8 still contains a considerable reservoir of molecular gas and the dense gas that 
has not yet been dissipated by UV radiation, it is likely to undergo star formation activities. 
The cloud NMA-8 can be at a stage of ongoing massive star formation as well as dense gas 
formation. 



4.2.2. NMA-9 

This is the second largest cloud (size of ~ 28 pc and mass of ~ 6.4 x 10^ M0) with the 
lowest SFE. Its low SFE24^m suggests a lack of recent massive star formation. This cloud 
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includes dense molecular gas, which is suggested by the detection of HCN line emission. 
Therefore, NMA-9 might be just forming dense molecular gas in the GMC but not be led to 
form stars yet. 



4-2.3. The other clouds 



The other clouds show much lower values in SFE24^m than NMA-8 and NMA-4. We are 
not sure whether they include dense gas or not because HCN emissions were not detected 
over Scrims in some clouds (NMA-1, NMA-6, NMA-7, and NMA-10). The S'hcn/5'co ratio 
upper limits a llow for the possibility that CO clouds h ave dense gas fraction similar to 
typical GMCs faelfer fc Blitzlll997bl : fBrouillet et al.lboosi ). Much deeper HCN observations 
are necessary to distinguish different properties from the dense gas fraction. Even though 
as for other clouds (NMA-2, NMA-3, and NMA-5), they have not been observed in HCN 
emission, we briefly describe the possibility of de nse gas formation i n the third largest cloud, 
NMA-2, since CO(3-2) emission was detected jTosaki et al.lboOTl l. The CO(3-2) is often 
used as a tracer of warmer and/or denser molecular gas [rin^ ~ lO^cm"'^, Tent = 33 K) than 
CO(l-O). In fact, CO(3-2) is a better tracer of density than temperature in the region where 
the CO( 3-2)/CO(l-0) ratio is less than 0.7 when assuming the kinetic temperature is low 
(~ 20 K; iMuraoka et aPboO?! ). Its CO ( 3-2) /CO (1-0 j rati o (~ 0.5) has a value similar to 
that of the second largest cloud NMA-9 (jTosaki et al.ll2007l ). Since NMA-9 shows no sign of 
star formation but dense gas formation in our results, the detection in CO (3-2) emission on 
this cloud could be due to the dense gas formation rather than heating by embedded young 
stars. Considering that NMA-2 has low SFEHa,24/im and shows little sign of star formation, 
NMA-2 is likely to have a similar status regarding dense gas and star formation, i.e., the 
clouds include a significant fraction of dense gas but have low star formation activities. 

In summary, NMA-9, NMA-8, and NMA-4, seem to have significant different properties 
in terms of dense gas formation and massive star formation. This variety can be explained as 
differences in their evolutionary stages, from dense gas formation to star formation. Cloud 
NMA-9 could be at the stage of "dense gas formation with no massive star formation" . Cloud 
NMA-8 could be at more evolved stage of "dense gas and ongoing massive star formation" 
than NMA-9. Cloud NMA-4 could be at the most evolved stage of "remnant" cloud where 
most of the molecular gas has already been converted into stars and been dispersed by newly 
born massive stars. 
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4.3. Sequential Star Formation in NGC 604 

As we mentioned in the previous section, we found a distinct difference in tlie pliysical 
states of star formation among tlie GMCs in the NGC 604 region and explained it as a 
difference in their evolutionary stages. Furthermore, these GMCs are well correlated with 
the spatial structure, and some of them are also with kinematics, of the Ha nebula including 
the expanding shells (Section l3.4p . The physical properties and the evolution of GMCs are 
straightforwardly expected to be affected by the Hll region excited by the stellar activities 
of OB clusters. The Shell B might have been forrn ed by the blowout through the molec- 



ular tun nel over Shell A (ITenorio-Tagle et al.l 120001 ) and is estimated to be younger than 



Shell A (lYang et al.lll996l ). which is good agreement with that NMA-8 is at a younger stage 
than NMA-4 because NMA-4 has already dissipated parent molecular cloud as mentioned 
in Section 4.2. These results are explain ed with the scenari o of sequential star formation in 



NGC 604 that was previously proposed flTosaki et al.l 120071 ) 



Figures [7] shows the radial distribution of L(Ha)corr5 M^^j, and SFEHQ,24/xm for each 
cloud as a function of projected distance from the center of the main cavity in NGC 604 
(0l'^34™32!5, +30°47'0'.'5), where the central clusters are crowded most. The gradual decrease 
of L(Ha)corr as the radial distance increases is found in Figure [7] (top), except the NMA- 
8 region (Shell B region). On the other hand, we found the gradual increase of MJj^^j, 
which implies that the star formation are active enough to consume and/or to dissipate 
parental clouds. We found a clear radial dependence in SFEHa,24^(m; the closer to the central 
clusters the GMCs are located, the higher SFE they have. We should note that this is 
because of the combination of the radial dependency of both L(Ha)corr and Mj^^j rather 
than only of L{B.a)coTi distribution. Besides, the total flux of L(Ha)corr is considered to be 
contributed from both the central clusters and the young stars that might be embedded in 
the clouds (Section I1X2D- As for NMA-2, NMA-6, and NMA-9 we don't know whether 
they are physically associated with Ha emission or not. In the case that they are observed 
at associations by chance, the values of SFE that we derived give just the upper limits. 
Nevertheless, the SFE on the bottom panel in Figure [7| still shows tendency that the SFE 
decrease as the distance from the central clusters increase. 

Figure [8] also shows the radial distribution of Sucn/Sqq and SFEggcHz as a function 
of projected distance from the center of the main cavity in NGC 604. We found a radial 
dependence in SFEggcHz, which are in good agreement with the behavior of SFEHQ-,24Aim (the 
bottom panel of Figure [7]) , even though these values were derived from different wavelength 
data. This implies that 89 GHz continuum can be another tracer of massive-star-forming 
region. On the other hand, the radial dependency of the dense gas fraction Suc^/Sqq is 
not clear from our data where we have estimated upper limits. The dense gas fraction 



- 18 - 



5'hcn/'S'co ill Figure [8] (bottom) is averaged over the extent of the CO clouds (denoted as 
circles). However, our high resolution images of CO and HCN emissions in Figure H] notify 
us that the hne ratio of Sucn/Sqq shows a gradient within a single cloud because two HCN 
components were detected at the northern part of the CO clouds NMA-8 and NMA-9 rather 
than the peak of the CO emission (see Section 13. 2p . The line ratios of 5'hcn/5'co (and 
its locations) that were measured over the region where HCN emissions were detected are 
represented as an open square in the bottom panel of Figure [HI Their values are higher than 
the averaged value which are represented as black circles in Figure M and their locations are 
close to the direction of the central clusters as compared with the line ratios within a single 
cloud. Thus from this figure it is found that the dense gas formation has been occurred at 
the near side of a cloud toward the central clusters. 

In the light of our high-resolution observational results, we reconsider the "sequential 
star form ation" scenario in NGC 604 here in some more detail than the previously proposed 



scenario (ITosaki et al.ll2007l ). First, central clusters were formed, whose powerful stellar wind 



m ade the main cavity with Shell A; this is referred to as "first-generation star formation" 



m iTosaki et al.l ( 120071 ). Inside Shell A, r nolecular gas around the ce ntral part of the stellar 
clusters is almost cleared out (Section 3.4: lMafz-Apellaniz et al. I l2004h . The "remnant" cloud 



i.e., NMA-4, is on the rim of Shell A. It is suggested that NMA-4 is a cloud swept up by the 
stellar wind from the central stellar cluster. The second-generation stars are formed in the 
swept-up cloud through dense gas formation and such new stars ionize the surfaces of the 
parental clouds. 

The cloud with "dense gas and ongoing massive star formation", i.e., NMA-8, is located 
on the edge of Shell B, while molecular gas inside Shell B is almost cleared out. Similar to 
the correlation between Shell A and NMA-4, NMA-8 could be formed with the ambient gas 
swept-up on Shell B and compressed by stellar winds. In such a compressed dense gas on 
Shell B, dense gas, and then stars are newly formed as is evidenced by that dense gas and 
star forming regions are located at the inner wall of NMA-8 toward the central stellar cluster. 
Besides, the star forming region in the NMA-8 is much closer to the central stellar cluster 
than dense gas forming region, and hence, central cluster, massive star forming region, dense 
gas forming region, and molecular cloud are aligned in this order. The newly born stars in 
NMA-8 associated with Shell B are the latter generation stars of sequential star formation 
which are induced by the prev ious generation stars; they are "second-generation stars" as 



referred to iTosaki et al.l ( 120071 ). or latter generation stars. 



On the other hand, the "dense gas forming" cloud, i.e., NMA-9, is located on the 
edge of the extended diffuse Ha emission up to a few hundreds of parsecs south of the 
central cluster. The dense gas forming region which is associated with NMA-9 cloud, are 
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located to the direction of central stellar cluster rather than the emission peak of NMA-9 
cloud. The expansion of the H II re gion into surrounding me dium can trigger star formation 



through dense gas formation (e.g., lElmegreen fc Lada]ll977l ). The NMA-2 is located in the 



north of the diffuse Ha nebula, while NMA-9 is in the south of nebula (Section 4.2). Their 
projected distances from the central stellar cluster are almost the same. If we suppose that 
the propagation initiated from the central stellar cluster is isotropic, dense gas is expected 
to be formed in NMA-2 as well as in NMA-9. The NMA-9 might be compressed by the 
expansion of the H II region. 

In summary, we explain the following scenario of sequential star formation in NGC 604. 
The first-generation stars (the central clusters) are formed in a GMC. The central clusters 
swept up the matter to form Shell A, where dense gas was formed by compression of the 
cloud material. Then second-generation stars form in the dense gas, which have already 
been dissipated in the parental clouds (the "remnant" cloud; NMA-4). Shell B is formed by 
blowout through the molecular tunnel, where another star formation occurs on the "on-going 
dense gas and star formation" cloud (NMA-8). At the far end beyond these shells dense gas 
is formed in the quiescent molecular clouds (NMA-9) by the effect of the expansion of the 
H II region. 



5. Summary 

We presented the highest sensitivity and the highest resolution ^^CO images of the 
giant H ll region NGC 604 in the nearby galaxy M 33. Moreover, we also obtained the HCN 
and 89 GHz continuum images. This is the first detection of HCN and 89 GHz continuum 
emission toward NGC 604 reported in the literature. 

1. We identified ten CO molecular clouds that have sizes of 5 to 29 pc and masses of (0.8 
- 7.4) xlO^M©. They have comparable sizes and masses to the typical GMCs in our 
Galaxy. 

2. The HCN-to-CO flux ratio of NMA-8 is 0.024 ± 0.007 and that of NMA-9 is 0.018 ± 
0.007. These values are comparable to those of typical GMCs in our Galaxy and in a 
neighboring galaxy M31. 

3. We derived the SFEh„. 24^m for each cloud from dust-corrected Ha luminosity (com- 
bination of Ha and 24 fim luminosity). We found that the SFEHQ,24/xm "was in good 
agreement with that from 89 GHz continuum. 
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4. Three CO clouds show a good spatial and velocity correlation with the bright Ha 
shells. The CO clouds show variations in SFEs, as estimated from the 89-GHz and 
combination of Ha and Spitzer 24 /xm luminosity. These variations of their properties 
could be explained as a snapshot of molecular clouds at different evolutionary stages. 

5. A clear radial dependence of SFEs from the center of the main cavity was found. 
Besides, we found that the HCN and 89 GHz continuum emission show offset from the 
CO peak and are distributed to the direction of the central cluster. We interpret our 
results as a scenario of sequential star formation in which the massive star formation 
propagates through the expansion of the Ha emission nebula excited by the central 
OB star cluster. 
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Fig. 1. — Positions of the five fields of view (FoVs) observed in CO (1-0) overlaid on the 
distribution of the Ha which is available in the HST archive. The archival HST WFPC2 
image of NGC 604 from program No. 5237 (PI J. Westphal) is used here. Each circle 
corresponds to the NMA primary beam (FWHM= 59") at 115 GHz. The field names are 
also labeled. A bar on the right bottom in the figure represents linear scale of 100 pc at the 
distance of M33. 
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Fig. 2. — Total integrated intensity map of ^^CO over five FoVs. Primary beam correction 
was applied to this image. Contour levels are -3, 3, 6, 9, 12, 15, 18, and 21crrms, where 
1 CTrms = 1.8 X 10^ M Jy sr~^ km s~^ . To convert to Jysr~^, the conversion factor 2.8 10"^*^ 
srbeam"^ was used. The crosses represent the positions of the 12 identified clouds and their 
names are also labeled. The circles represent five FoVs that we observed in CO line emission. 
The parts of the emission outside the FoVs are masked out. 
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Fig. 3. — The spectra of the two HCN components. The left panel is one component 
associated with the CO molecular cloud NMA-8, while the right panel is the other associated 
with NMA-9. The CO spectra for NMA-8 and NMA-9 are also shown, which are measured 
in the same region on the convolved image with the same synthesized beam and with the 
same velocity resolution as the HCN image. 
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Fig. 4. — Total integrated intensity map of HCN (red contours) and CO (black contour) 
overlaid on the Ha image (gray scale) for NMA-8 (left) and NMA-9 (right). The inte- 
grated velocity ranges are from Vlsr = — 245kms~^ to — 252kms~^ for NMA-8, and from 
Vlsr = — 222kms~^ to — 225kms~^ for NMA-9. The contour levels for HCN map are 
— 2.5, 2.5, 5 (Tj-nis) where Icrj-ms — 1-5 x lO^MJysr^^kms"^ The conversion factor from 
Jybeam"^ to Jysr~^ is 9.3 x 10^^*^ srbeam"^ for the 6'.'6 x 5'.'3 resolution data. The negative 
contours are drawn as gray contour. The black contour levels for CO are 3, 6, 9, 12, 15, 
18, and 21 arms- The gray circles with the diameter of 77" are two FoVs that we observed 
in HCN and 89 GHz continuum emission. The black arrows indicate the region where HCN 
emission was detected over two channels and with > 3 a. The parts of the emission outside 
the FoVs are masked out. Primary beam correction was applied to these images. 



- 28 - 




Fig. 5. — (ieft) The 89 GHz continuum map (red contour) and CO (black contour) overlaid 
on the Ha (gray scale). The red contour levels for 89 GHz continuum map are -2, 2, 4, 
and 6(Trms, where 1 a^ms = 0.6MJysr"^, corresponding to 0.6 mJybeam"^ for the 6'.'9 x 5'.'5 
resolution data. The negative contours are drawn as gray contour. The black contour levels 
for CO are 3, 6, 9, 12, 15, 18, and 21 aj-ms- The gray circles with the diameter of 77" are two 
FoVs that we observed in HCN and 89 GHz continuum emission. The parts of the emission 
outside the FoVs are masked out. Primary beam correction was applied to this image, (right) 
Comparison of 89 GHz continuum emission (red contour) with 24 /im emission (gray scale) 
overlaid CO image (black contour). The 89 GHz source names (I-IV) are labeled. The cross 
symbol indicates the position of the center of the main cavity in NGC 604. The green star 
symbo ls represent the peak position of the 8.4 GHz components (A-F) by lChurchwell fc Goss 
f ll999h . 
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Fig. 6. — (left) Interferometric integrated intensity maps of ^^CO(J=1-0) emission (cyan), 
HCN emission (red) and 89 GHz continuum (yellow) overlaid color image of Ha (gray scale) 
by HST, which are the same images as Figure 1-5. But yellow contour levels are > 4 a^-ms and 
red contour levels are 4, 5 (Jrms in order to show clearly the peak positions. The gray contour 
levels are 0.2 % of the peak intensity so as to show the extent of the diffuse Ha emission 
nebula. The pink star symbol indicates the position of the center of the main cavity in 
NGC 604. The gray circles with the diameter of 77" are two FoVs that we observed in HCN 
and 89 GHz continuum emission. The cloud numbers and linear scale bar are also plotted. 
(right) The trimmed image of Figure 1 draws the features of mains shells and annotations 
discussed in the text. The locations of two shells named A and B are indicated on the map. 
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Fig. 7.— L(Ha)corr (top), M*^^^ {middle), and SFEhc, 24^m {top) distribution as a function 
of projected distance from the center of the main cavity in NGC 604. Its coordinate of 
(0l'^34™32!5, +30°47'0'.'5) is indicated as a pink star symbol in Figure El The cloud numbers 
are also plotted. In the top panel, the location of two Ha shells are indicated for reference. 
Note that common errors are excluded in order to compare each clouds relatively. 
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Fig. 8. — SFEgQGHz (top), and Shq^^/Sqq {bottom) distribution as a function of projected 
distance from the center of the main cavity in NGC 604. Its coordinate of (0l''34'^32!5, 
+30°47'0'.'5) is indicated as a pink star symbol in Figure El The filled circles represent the 
detected clouds, while the open circles show the upper limit for the clouds where were not 
detected in HCN and 89 GHz continuum emission. The open squares represent the line ratio 
that measured in the region where HCN emissions were detected. The distance of these two 
plots were measured at the peak of HCN emission. The cloud numbers are also plotted. 
Note that common errors are excluded in order to compare each clouds relatively. 
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Table 1. ^^CO(J=1-0) Observational Properties 



FoV 


Field Center 


Array 


Integration 


Beam Size 


P.A. 


^ rms 


Flux Ratiol^ 




(RA.DEC) 


Configuration 


Time (hours) 


(arcscc) 


(degree) 


(MJy sr-i) 




FOVl 


(0l'^34'°31?7, +30°47'41'.'0) 


C,D 


8.2 


5.4x4.5 


-38.1 


75 


0.51±0.09 


FOV2 


(0l''34™32!7, +30°47'13'.'9) 


C,D 


7.8 


5.4x4.5 


-38.1 


75 


0.23±0.04 


FOV3 


(0l''34™33';3, +30°46'47'.'4) 


C,D 


6.0 


4.1x2.6 


-46.7 


147 


0.62±0.11 


FOV4 


(0l''34™30?9, +30°46'53'.'7) 


C,D 


7.9 


5.4x4.5 


-38.1 


75 


0.32±0.06 


FOV5 


(0l''34™34!5, +30°46'18'.'4) 


C,D 


6.6 


4.1x2.6 


-46.7 


147 


0.63±0.11 



''The rms noise was measured in the emission-free region with 1-MHz resolution on the primary-beam-corrected images. The 
unit was converted from Jybeam"^ to Jy sr~^ using ^hea,m = 4 hi 2 ^™^maj sin^niin 1 where fiboam is the synthesized beam in 
steradian, ^maj a^nd ^min ^^re the synthesized beam in arcsec. The conversion factors arc 6.5x10"^'' srbcam"^ for the 5'.'4 X 4'.'5 
resolution data, and 2.8x10"^*^ srbeam"^ for the 4'.'1 X 2'.'6 resolution data. 

'^The ratio of interferomctric fiux to the single dish flux (R. Miura et al., 2010, in preparation), which is measured in the box 
region with 42"x42"at the center of each FoV. 



Table 2. Observational parameters of three probes 



Parameter 


12CO(1-0) 


HCN(l-O) 


89 GHz Continuum 


Frequency (GHz) 


115.271202 


88.6304157 


88.9790432 ±0.512 (LSB)'' 


Observing Period 


2004 Nov - 2006 Doc 


2004 Nov - 2006 Jan 


2004 Nov - 2006 Jan 


Observed Field 


FOVl, 2, 3, 4, 5 


FOV3, FOV5 


FOV3, FOV5 


Array Configurations 


C, D 


C, D 


C, D 


Correlator 


FX 


FX 


UWBC 


Bandwidth (MHz) 


32 


32 


1024 


FWHM of FoV (") 


59 


77 


76 


Beam Size (") 


(see Table 1) 


6.6x5.3 


6.9x5.5 


Velocity Resolution (km s^-*-) 


2.6 


3.4 




Integration Time (hour) 


(see Table 1) 


19 


19 


o-rms (MJy sr-i)*^ 


(see Table 1) 


21 


0.6 



''The frequency range corresponding to the line emission as HCN and HCO+ (89. 188526 GHz) were removed 
from the original visibility data. 

''The conversion factors are 9.3x 10^'" sr beam~^ for the 6'.'6 x 5'.' 3 resolution data, and l.Ox 10^'-*sr beam"'^ 
for the 6'.' 9 X 5" 5 resolution data. 
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Table 3. Cloud Properties for Clouds Found in M 33 



Cloud 


Alternative'' 


a(J2000), <5(J2000)'' 


Vlsr" 


Vfwhm'* 


Sco " 








No. 


Name 




(kms-i) 


(kms-i) 


(Jy kms""*^) 


(pc) 


(105 A/0) 


(lO^Mo) 


NMA-1 




0l'"34™31!59, +30°46'54y88 


-241 


4.9 


14.4±1.5 


< 9.4* 


1.6 


< 0.3 


NMA-2 


EPRB-101 


0l'^34™31?59, +30°47'42y03 


-244 


10.4 


22.3±1.6 


19.2 


2.5 


2.9 


NMA-3 




0l'"34™31!94, +30°47'25yi4 


-250 


7.3 


6.6±1.0 


< 4.9* 


0.8 


< 0.4 


NMA-4 




0l''34™32!29, +30°46'57y40 


-247 


8.4 


6.8±1.0 


10.1 


0.8 


1.0 


NMA-5 




0l'^34™32?74, +30°47'39yi8 


-242 


7.0 


10.7±1.7 


< 10.9* 


1.2 


< 0.9 


NMA-6 




0lh34™32;79, +30°46'13y55 


-240 


7.9 


18.4±2.2 


21.9 


2.1 


1.8 


NMA-7 


WS-1 


0lh34™33?25, +30°46'30yi9 


-249 


7.6 


18.6±1.7 


19.2 


2.1 


1.5 


NMA-8 


WS-2, EPRB-8 


0l'^34™33?42, +30°46'47y03 


-243 


10.2 


65.2±2.7 


28.9 


7.4 


4.1 


NMA-9 


WS-4, EPRB-9 


0l'^34™34!47, +30°46'20y00 


-221 


9.6 


56.7±2.4 


27.9 


6.4 


3.6 


NMA-10 




0l'^34™35!16, +30°47'0iy50 


-228 


4.8 


27.2±3.2 


11.0 


3.1 


0.4 


NMA-11 




0l''34™34!07, +30°47'09y47 


-239 


<2.6 


12.4±1.4 


< 20 


1.4 


<0.7 


NMA-12 




0l''34™34?64, +30°46'35yi0 


-228 


<2.6 


10.7±1.2 


< 20 


1.2 


<0.7 



Note. — The derived properties of identified clouds were described shortly in text. Approximate errors are (5VpwHM = 1.3 km s ^ 
(i.e., one-half the velocity resolution), SM^^l = 38%, and SM^ir = 40%. 

'^WS-n: from lWilson fc Scovillj l|l992h . EPRB-n: from lEngargiola et alj \200ch . 

''The positions of clouds wore obtained from the position of the peak of emission. 

'^The positions along the velocity axis were obtained from the position of the peak of the observed (Gaussian-fit) LSR velocities. 
'^Thc full width at half-maximum velocities are deconvolved with a velocity resolution of <5Vfwhm = 1.3kms~*. 
"The flux densities were corrected primary beam attenuation. 

*D = Da + Dg, where Da and Dg are the deconvolved full width at half-maximum diameters in right ascension and declination 
directions. Typical error of the cloud size differs depending on where the clouds are located because the synthesized beam size differs on 
each FoV. Typical errors are 5D = 5.1pc for FOV3 and FOV5, 5D = 6.6 pc for FOVl, FOV2 and FOV4. 

^Since the cloud size in right ascension of NMA-1, NMA-3, and NMA-5 arc less than beam size in the same direction, the cloud size 
D is given as the upper limit. 
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Table 4. Cloud Properties and Star Formation Activities for Clouds 



Cloud 






Shcn/S'^q'^ 


SggGHz'^ 


SFEgsGHz'^ 


logLcorr(J?a)° 


SFEHc«,24Mm'' 


No. 


(pc) 


(Jy km s^"*") 


(10-2) 


(mjy) 


(%) 


(ergs-1) 


(%) 


NMA-1 


86 






<2.6 


<5.0 


38.58 


6.2±3.3 


NMA-2 


214 










37.53 


0.4±0.2 


NMA-3 


120 










37.78 


2.5±1.3 


NMA-4 


35 






3.4±0.9 


12.4±6.6 


38.75 


17.2±9.3 


NMA-5 


152 










37.64 


1.1±0.6 


NMA-6 


246 


<0.3 


<2.8 


<2.1 


<3.1 


37.32 


0.3±0.2 


NMA-7 


194 


<0.3 


<2.7 


<1.8 


<3.6 


37.96 


1.6±0.8 


NMA-8 


138 


1.1±0.3 


2.4±0.7 


9.1±1.5 


5.1±2.3 


39.08 


5.9±3.0 


NMA-9 


320 


0.6±0.2 


1.8±0.7 


<1.8 


<1.2 


37.67 


0.3±0.1 


NMA-10 


261 


<0.4 


<2.8 


<2.3 


<2.3 


38.08 


1.1±0.6 


NMA-1,4 


60 


<0.4 


<3.2 











''The projected distance from the central cluster, dc, is given in parsec relative to the position of 
(0l'>34™32';4, +30°47'00('3) 

''For the clouds where no HCN emission was detected, the upper limits of line ratio are given by using 
the detection limit of 3 cTrms . 

'^The CO flux densities, Sqq, were measured over the extent of CO emissions (> 3 a-ims) on the convolved 
images with HCN synthesized beam. 

''The SFEgg GHz was obtained using the following equation, SFEgg gHz = Afstar/(Afstar + A/j^Qj), where 
-^mol represents the missing-flux-corrected molecular masses. See details in the text. 

''The extinction-corrected Ha luminosity is defined as a line ar combination of th ese two luminosities: 
icorr(Ha) = Lobs(Ho)-f axL(24^m), where a = (0.031±0.006) jCalzetti et al.ll2007l) . 5Lcoo (Ha) ~ 28%. 

''The SFEHc<,24fjm wS'S obtained using the same equation as SFEgg ghz- See details in the text. 
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Table 5. Star Formation Efficiencies of nearby Hll regions 



Cloud 


Mass (lO'^ M0) 


Q* (lo'^-'s-i; 


) SFE 


References 


Orion A 


1 


2.7=" 


2.2 


1 


W49A 


7 


182 


18.2 


1, 2, 3 


30 Doradus 


9b 


1122<: 


52.8 


2, 4 


NGC604 


51^ 


478° 


7.4 


this worlc 


...NMA-4 


2d 


28 


12.5 


this work 


...NMA-8 


12^ 


75 


5.1 


this worlc 



Reference s. — ( llWimams fc McKej lll997^ and references 
th erein; (2lKennicuttl lll984l' l: (3 lWard-Thompson fc RobsonI \l99(t) : 
(4 lFukui et ahl l l200Sf) 

Note. — The SFE was calculated assuming the same IMF (2.35) 
throughout all regions. The typical error of SFE is ~50%. 

"The ionizing photon rates in O rion A is determined from the stellar 
types in lCardeUi fc ClavtonI lll988l '). 

'^This value is summation of mo lecular cloud mass of GMCs which are 
annotated by 183, 186, and 189 in lFukui et al.l ||2008|) . which are located 
withi n the extent of Ha emission in 30 Doradus (~370 pc; iKennicuttI 
Il984l '). 

'^The ionizing photon rate Q* is derived only from Ha luminosities 
l lKennicuttlll984h , and hence the SFE can give lower limit. 

'^They represent the missing-flux-corrccted molecular masses (M^^^). 

"The ionizing photon rate Q* is derived from the extinction-corrected 
Hq luminosities measured over a whole extent of Ha. 
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